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Abstract. - Aging in spin glasses is analyzed via the Probability Density Function (PDF) 
of the heat transfer over a short time St between system and heat bath. The PDF contains 
a Gaussian part, describing reversible fluctuations, and an exponential 'intermittent' tail. We 
find that the relative weight of these two parts depends, for fixed St on the ratio of the total 
sampling time t to the age t w . Fixing this ratio, the intensity of the intermittent events is 
proportional to 5t/t w and independent of the temperature. The Gaussian part has a variance 
with the same temperature dependence as the variance of the equilibrium energy in a system 
with an exponential density of states. All these observations are explained assuming that for 
any t w , intermittent events are triggered by local energy fluctuations exceeding those previously 
occurred. 



Introduction. Aging glassy systems undergo reversible equilibrium-like fluctuations 
alongside with irreversible configuration changes. As shown experimentally in [1-3] for a 
number of different cases, these two aspects can be disentangled by a statistical analysis of 
a mesoscopic noisy signal: within the probability density function (PDF) of the signal, the 
fluctuations give rise to a Gaussian peak, and the irreversible changes to an asymmetric 'inter- 
mittent' tail. This probe of aging dynamics considerably differs from approaches utilizing the 
fluctuation-dissipation theorem, and raises new theoretical issues regarding e.g. the statistical 
characterization of the intermittent events and their relationship to the thermal fluctuations. 
We consider the statistics of the heat-flow over a small interval St. As was the case in [1-3], 
two sub-processes are discernible in the PDF of the data: (i) intermittent events, referred 
to as quakes, which appear in the left tail of the PDF. These irreversibly release the excess 
energy trapped by the initial quench, (ii) Gaussian fluctuations which occur between the 
jumps of type (i) and correspond to reversible energy fluctuations within thermalized attrac- 
tors. Type (i) events are rare compared to type (ii), but are nevertheless orders of magnitude 
more frequent than Gaussian fluctuations of similar size. Importantly, the rate of intermittent 
events decreases with the reciprocal of the age t w and is independent of the temperature T. 
Conversely, the Gaussian fluctuations are t w independent, but strongly T dependent. 
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Fig. 1 - The PDF's (dots) of the heat transfer H(8t,t,t w ), sampled in the observation window 
[t w ,t w +t] = [40000,42000] for St = 10,20,40,60 and 80. The temperature is T = 0.3. Insert: the 
mean (circles), variance (diamonds) and skewness (stars), are plotted versus the time difference St. 

Fig. 2 - Data collected as in Fig. but with the observation window [t w ,t w + t] = [2000,4000]. A 
strong asymmetric tail appears, producing a large negative skewness when the observation time t and 
the age t m are of comparable magnitude. 



In this Letter we identify the intermittent events with the large and irreversible configu- 
ration changes of glassy systems already described in ref . [4] , where they were similarly called 
quakes. We hence borrow the idea that their occurrence is due to energy fluctuations larger 
than any hitherto occurred, i.e. fluctuation records and show that record-induced dynam- 
ics [4,5] predicts the key dynamical features of the heat transfer PDF. The predictions are 
first spelled out and tested against simulation data for a prototype glassy system, a three di- 
mensional Ising spin glass with short-ranged couplings. We then argue in greater detail for the 
theoretical link between record-sized energy fluctuations, density of states within metastable 
attractors and intermittency, specifically addressing EA spin glasses and other models with 
short ranged interactions. We conclude with a summary and a brief discussion. 

Measurable consequences of record-induced dynamics. Non-equilibrium aging is de- 
scribed in Ref. [4] in an cartoon- like fashion as a series of irreversible 'quakes' between dy- 
namically inequivalent metastable attractors of an energy landscape. The underlining physical 
assumptions are: Firstly, that the stability of the 'current attractor', as gauged by an 'exit' 
barrier, increases by a tiny amount with each quake. Hence a positive energy fluctuation larger 
than all previous ones, i.e. a record, triggers the irreversible change of attractor. Secondly, 
that de-facto irreversibility of the quakes arises due to the large energy loss. I.e., crossing 
extremal barriers triggers irreversible configuration changes, while smaller barriers are crossed 
reversibly. Thirdly, saddles are visited rarely and in a statistically independent fashion. With 
these assumptions, the quakes occurring in the interval (t w ,t w + t] are well described by a 
Poisson process with logarithmic time arguments [4,8], and, by our identification, so are the 
intermittent events. 

The average number of intermittent events is thus 

(ni(t w ,t w +t)) = cdog(l +t/t w ), a>0, (1) 
where the parameter a slowly increases with the system size but is independent of T [8]. 



Paolo Sibani and Henrik Jeldtoft Jensen: Intermittency and aging 



3 



Two predictions which are measurable from the PDF of intermittent data without access to 
microscopic information, are easily derived from Eq. ^ 

Firstly, for St <C t w , a single event occurs in the interval St with probability 

P{tw , tw+ 5t](l) = a^+0(5t/t w ) 2 , (2) 

while multiple events have zero probability to linear order in 5t/t w . Thus, the rate of inter- 
mittent events decays as a/t w and is independent of the temperature. 

Secondly, the number of intermittent events falling in an arbitrary time interval t grows 
on average with log(l + t/t w ). As reversible energy fluctuation occur at a constant rate, their 
number uq in the same interval grows linearly in t. The ratio nj/riQ is hence given by 

— oclog(l +t/t w ). (3) 
n G t w 

The statistical properties of records in a sequence of random numbers are largely insensitive 
to the distribution from which these numbers are independently drawn [4,5]. This leads to the 
predicted T independence of the intermittency rate, which is however a puzzling conclusion 
when thermally activated barrier climbing arguably produces the attractor changes. In the 
penultimate section of the Letter we analyze this issue in some detail, define the attractors 
more concretely for a system with short-ranged interactions, and derive Eq. El The main 
result is that the cancellation of strong T dependences needed to ensure consistency between 
record-dynamics and thermal barrier climbing occurs precisely if the energy fluctuations are 
drawn from an exponential distribution. In our case, this distribution coincides with the 'local 
density of states' (LDOS) characterizing the energies available in local thermal equilibrium. 
Writing the LDOS as 

V(e) oc exp(e/e ) (4) 

produces the heat capacity 

Ct ,(T)cxe2/(e -T) 2 , (5) 
and leads to a Gaussian variance of the heat transfer data 

a% cx c v (T)T 2 cx f^V) , T < e . (6) 

This above equation is our third measurable prediction. 

Simulation results. The heat-flow data stem from simulations of a standard 16 3 
Edwards- Anderson spin-glass with nearest neighbor Gaussian couplings [6]. For speed we 
employ an event driven simulation technique [7,8], whose 'intrinsic' time corresponds, for 
large systems, to the number of Monte Carlo sweeps. The Boltzmann constant is set to one, 
and age zero is fixed by the instantaneous deep thermal quench preceding the simulations. 

We denote by E{t w ) the total energy at age t w , pick St < t w and consider a series of n 
contiguous time intervals of length St. For the k + l'th interval, k = 0, 1 . . . n — 1, the heat 
H exchanged is calculated as the difference E(t w + (k + l)St) — E(t w + kSt). Accordingly, H 
is negative when heat leaves the system. Sampling H over a total observation time t = nSt 
produces the raw data behind the Probability Density Functions (PDFs) discussed below. 
Since, as it turns out, the form of the PDFs depends on the three parameters St, t and t w , the 
heat transfer is denoted H(5t,t,t w ). Omitting the parametric time dependences, the PDF of 
H is denoted by pu{x), whereas capital letters e.g. P q are used for cumulative distributions. 
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Fig. 3 - The probability P q (e) that the amount q of energy intermittently released is larger than 
e, is plotted versus e for ages t w = 1000, 2000, 4000, . . . 64000. Lower insert: (unnormalized) PDF 
of the heat transfer H(0.01t w , t w /2, t w ). The full lines — indistinguishable on the scale of the plot — 
are fits to Gaussians with zero average. Upper insert: the relative weight ui/ug (see text) of the 
intermittent tail is obtained varying t for fixed St from seven values of H(0.005t w , t, t w ) with t w = 4000 
and t = 200, 600, 1000, 2000, 3000, 4000 and 6000. The straight line is predicted by the theory. The 
simulation temperature is T = 0.3. 

Fig. 4 - Left frame, symbols: reciprocal slope of — log(P 9 (e)) for large e, see Fig. [3] as a function of 
T. The dotted band is proportional to T -0 ' 2 . Right frame, symbols: variance a G of the Gaussian 
part of the PDF. Full line: theoretical prediction a G oc T 2 c v (T), with c v (T) given by Eq.0 



To improve the statistics simulations are repeated for 625 or, mainly, 1250 independent runs 
for each set of physical parameters. 

Figures □ and H show data taken for t = 2000, T = 0.3 and five different 5t's. In Fig[T] 
t w = 20t while t w — t in Fig. O In agreement with previous observations [1-3,9], the 
PDF's have a Gaussian part and an exponential tail. In our figures, this tail becomes more 
pronounced as t t tw, a characteristic shift from equilibrium to non-equilibrium dynamics 
usually studied through correlations and response functions [10-12]. 

The asymmetry of the tail shows that large negative H values are associated with irre- 
versible energy losses. These events alone carry the linear St dependence of the mean, variance 
and skewness shown in the inserts of Figs. and [21 i- e - the Gaussian part of the PDF has 
always zero average and a variance which depends on T but not on St or t w , see e.g. Fig. [2 
(NB: the scale differs from FigQ] ). The linear growth of the variance, also shown in the 
inserts, continues up to St t w (not shown) and indicates statistical independence of the heat 
exchange over consecutive time intervals. 

For fixed t/t w and increasing t w , the rate of intermittent activity decreases as l/t w in a 
temperature independent fashion. This l/t w scaling is shown for T = 0.3 in Fig. [31 where a 
64-fold increase of t w is offset by choosing St = t w So with <5o = 1/100. The seven t w values 
considered are 1000 • 2™, for n = 0, 1, . . . 6. 

The good agreement with the prediction given in Eq.[21is demonstrated in the upper insert 
of Fig.[3by plotting (nj /na)(t/t w ) versus log(l +t/t w ). For fixed St — 20 and t w = 4000, we 
consider observation times t = 200,600,1000,2000,3000,4000 and 10000 (circles). The line 
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is a linear fit to the data. The unnormalized PDF's of H(t w /100,t w /2,t w ) (dots) and the 
respective Gaussian fits (overlapping lines) are plotted in the lower insert. The main panel 
shows the 'tail' distributions P q (e) . For large negative H values these describe the difference q 
between the empirical PDF value and the Gaussian probability for an energy fluctuation of size 
H(t w So, t w /2, t w ). The latter is estimated by fitting a Gaussian of zero average to the 'central 
part' of the PDF, as delimited by a symmetric cutoff at w ±2<tg- Pq( e ) is later interpreted as 
the probability that a single quake releases an energy q larger than e. The size distribution 
of the quakes is nearly exponential, with a decay constant l/qo- The St/t w scaling procedure 
yielding Fig. 02 was repeated for T = 0.2,0.4 and 0.5, with the same result. Additionally, for 
T = 0.25, 0.35, 0.45, six sets of simulations were carried out at fixed t w = 16000. For each T, 
averaging over the outcomes yields the qo and Oq values plotted in Fig. 0] as circles. The full 
line in the right panel of Fig. 01 is obtained by Eqs.ElandEl using a fitted pre-factor, which 
amounts to a vertical parallel shift of the theoretical line. The value eo = 0.86 used in Eq. [S] 
is very close to reported values of the critical temperature in 3d Gaussian spin glasses [13]. 

We note that while the statistical variation of Oq over different simulations set is insignif- 
icant, qo is estimated with decreasing accuracy as T and ac increase: the distinction between 
quakes and fluctuations is moot for e smaller than ps 2<tg, and the undetermined part of the 
distribution must be excluded from the fit. With the present numerical effort and accuracy, 
we can only analyze data below T = 0.6. The choice of cut-off gives a systematic error, not 
included in the ±<r error bars. The gray band shows the qualitative agreement with the T~ - 2 
scaling of the energy difference between local minima found in Ref. [8]. This temperature 
dependence is rather weak compared to the T dependence of Oq . 

Record fluctuations and intermittency. - In extended systems with short ranged correla- 
tions, thermal fluctuations are localized within well separated sub-volumes whose linear size 
is of the order of the thermal correlation length [14, 15]. As e.g. suggested in ref [16], each 
thcrmalized sub-volume can be treated as a small glassy system, with a 'local' energy land- 
scape. Such landscapes have numerous energy minima. For each minimum x, a metastable 
attractor can be defined as a configuration space neighborhood, e.g. the set of configurations 
connected to x by paths requiring energy changes below a given threshold [17]. The 'local 
density of states' is the density of state (DOS) restricted to the configurations belonging to the 
attractor. LDOS for several glassy systems [17] are indeed well approximated by exponential 
functions of the energy e measured from the lowest minimum of the attractor. The divergence 
of the heat capacity seen in Eq. [S] at T — eo implies that for T > eo the thermal stability 
of the attractor vanishes. In practice a rounding off occurs once the weight of configurations 
outside the attractor becomes important. 

The Arrhenius time belonging to the typical extremal barrier crossed at time scale t w is ob- 
tained as follows: first we note that the Boltzmann distribution corresponding to an exponen- 
tial LDOS is itself exponential, and given by Pe(6) = exp(— ba(T)), where a(T) = (1/T— 1/eo) 
is an effective inverse temperature. Secondly, with a constant microscopic attempt rate, the 
largest fluctuation on a time scale t w is the largest of 0(t w ) values drawn independently from 
Pe- Applying a well known mathematical result [18], b max (T,t w ) oc ]n{t w )/a(T). Hence, 
a(T) cancels, leading to the T independent Arrhenius time ta — exp(a(T)b max ) — t w . A 
small negative temperature shift imposed at t w generates a mismatch between the extremal 
barrier previously established and the thermal fluctuations available to cross it. The Arrhe- 
nius time and the intermittency rate acquire then a peculiar temperature dependence which 
has recently been verified in simulations [19]. 

Clearly, the sum of many independent energy fluctuations is a Gaussian quantity^), and so 



( 1 )Using a very small St the assumption of independence breaks down, and the Gaussian is replaced by a 
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is the central part of the heat transfer PDF, which arises from differences between equilibrium 
energy fluctuations. The Gaussian variance cr^ is therefore twice the equilibrium variance of 
the energy. The latter follows from Eq. 0] and textbook arguments of equilibrium statistical 
mechanics. 

Finally, figs lll3l show how the record dynamics scenario accounts for the dependence of 
the heat transfer PDF, on 6t,t and t w . Writing the latter as a a sum of a Gaussian and an 
intermittent part, pn ocpg +pi, the tail has the form 

pi{x) oc log(l + t/t w )^-p q {-x). (7) 

The factor p- combines the linear St dependence of the first three moments of H seen in 
FigsQ] and with the St/t w scaling of the intermittency PDF data shown in Fig.|3 for fixed 
t. The interpretation of p q (—x) as the PDF of the heat transferred in a single quake hinges 
on St/t w being proportional to the probability of a single intermittent event, see Eq|3 

An exponential p q could arise quite simply: the sequence of mainly downhill moves consti- 
tuting the quake releases an exponentially distributed energy if each move can be the last with 
equal probability. Lacking a prediction for the weak T dependence of the scale parameter go, 
we just note its qualitative agreement with the average energy difference between deep minima 
of the spin-glass landscape [8] . 

Discussion. The correspondence between intermittent events and attractor changes 
was previously investigated in ref. [9], using thermal quenches into inherent structures. We 
similarly link intermittency to quakes, and derive and verify novel analytical predictions for 
the heat transfer PDF, based on the assumption that record-sized positive energy fluctuations 
within the attractors induce the quakes. Record-induced dynamics has generic implications 
e.g.: the 6t/t w scaling form of the PDF exemplifies the decelerating dynamics found besides 
spin-glasses [11,12,15,20] in glasses [21], driven dissipative systems [22,23], and evolution 
models [24,25]. Nearly exponential densities of states are widespread [26,27], and were used 
already two decades ago [28] in mesoscopic models of glasses. We thus expect our analysis 
to broadly apply to glassy systems whenever reversible fluctuations and irreversible quakes 
are clearly separated, i.e. at sufficiently low temperature. With this limitation, it should also 
be extensible to other types of measurements, e.g. the fluctuations of two-point correlation 
functions [1] . Generalizations are needed to account for small but systematic deviations from 
t/t w scaling, i.e. sub-aging, whose origin is currently debated [29]. 

Extremal (Gumbel) distributions are ubiquitous, e.g. Chamon et al. [30] finds them to 
describe local correlators in aging systems, including systems without quenched disorder, and 
without local energy minima. Since the heat absorbed by an extremal equilibrium fluctuation 
is masked by the much larger amount of heat released in the following quake, the distribution 
of extremal energy fluctuations is not directly accessible through intermittent heat-transfer 
data. The relation with effective temperature descriptions of fluctuation-dissipation data [31, 
32] warrants a separate discussion. Here we only note that these descriptions concern the 
asymptotic limit of the fluctuation-dissipation ratio for large times. This regime differs from 
the one presently considered and is not always reached experimentally [2,12] 
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back-to back double exponential. 
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